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ABSTRACT 

Aims. We present a multi-wavelength analysis of the infrared dust bubble S24, and the extended IR sources G341.220-0.213 and 
G341.217-0.237 located in its environs, with the aim of investigating the characteristics of the molecular gas and the interstellar dust 
linked to them, and analyzing the evolutionary state of the young stellar objects identified there and its relation to S 24 and the IR 
sources. 

Methods. Using the APEX telescope, we mapped the molecular emission in the CO(2-l), *^CO(2-l), C'*0(2-l), and '^CO(3-2) lines 
in a region of about 5' x 5' in size around the bubble. The cold dust distribution was analyzed using submillimeter continuum images 
from ATLASGAL and Herschel. Complementary IR and radio data at different wavelengths were used to complete the study of the 
interstellar medium in the region. 

Results. The molecular gas distribution shows that gas linked to the S 24 bubble, G341.220-0.213, and G341.217-0.237 has velocities 
between ^8.0 km s“' and ^0.0 km s“', compatible with a kinematical distance of 3.7 kpc generally adopted for the region. The gas 
distribution reveals a shell-like molecular structure of ~0.8 pc in radius bordering the S 24 bubble. A cold dust counterpart of the shell 
is detected in the LABOCA and Herschel-SPIRE images. The presence of weak extended emission at 24 pm from warm dust and 
radio continuum emission projected inside the bubble indicates the existence of exciting sources and that the bubble is a compact Hit 
region. Part of the molecular gas bordering the S24 Hit region coincides with the extended infrared dust cloud SDC341.194-0.221. 
A molecular and cold dust clump is present at the interface between the S24 Hu region and G341.217-0.237, shaping the eastern 
border of the IR bubble. As regards G341.220-0.213, the presence of an arc-like molecular structure encircling the northern and 
eastern sections of this IR source indicates that G341.220-0.213 is interacting with the molecular gas. The analysis of the available 
IR point source catalogs reveals the existence of young stellar object (YSO) candidates linked to the IR extended sources, thus 
confirming their nature as active star-forming regions. Gas and dust masses were estimated for the different features. The 
totaf gas mass in the region and the H 2 ambient density amount to 10300 Mq and 5900 cm indicating that G341.220-0.213, 
G341.217-0.237, and the S 24 Hii region are evoiving in a high density medium. A triggering star formation scenario for the 
Hii region is investigated. 

Key words. ISM:molecules - stars:protostars - (/SAf. jlR dust bubbles - ISM:individual objects: S24 - ISM:individual ob- 
jects:IRAS 16487-4423 


1. Introduction 


IChurchwell et al.l(l2006ll2()()^ identified more than 600 infrared 
(IR) dust bubbles in the Spitzer 8.0 /tm-images of the Galactic 
Legacy Infrared Mid- Plane Survey Extraordinaire (GLIMPSE, 
iBeniamin et al.|[2()()^ between longitudes from -60° to h- 60°. 
These infrared (IR) dust bubbles are bright at 8.0 pm, show¬ 
ing the emission of polycyclic aromatic hydrocarbons (PAHs) 
typical of photodissociation regions (PDRs), and enclose 24 
pm emission, indic ating the existence of warm dust inside 
(I Watson et al.l l2008l) . These characteristics indicate the pres¬ 
ence of excitation sources (O- and/or B-type stars) that disso¬ 
ciate the molecular gas in the border of the PDRs and heat 
the dust through their UV photons. The dense shells surround¬ 
ing many IR dust bubbles are believed to favour the forma¬ 
tion of new generations of stars through the Radiatively Driven 
Implosion (RDI) or the collect and collapse (CC) processes 


(Pomares^t_alj2()()94^ peharveng et al.l200^ISamal et al.l2014t 

iDeharveng et al.ll2012ah . Indeed, detailed studies of these bub¬ 
bles have shown the presence of young stellar objects in 


their environs, although in some cases t riggered star formation 
could not be proved (se e for example [Alexander et al.l l2013t 
iDewangan & Oihall2013[) . 

We s elected the complex of IR dust bubbles S 21-S 24 iden¬ 
tified by IChurchwell et al.l (l2006h in the GLIMPSE images to 
perform a multiwavelength study of the interaction with their 
environs and search for young stellar objects in their vicin¬ 
ity. Here we report on the study of S 24. This bubble is cen¬ 
tered at RA,Dec.(J2000) = (16*52'"20', -44°28'10") or (l,b) = 
(341°203, -0°23) and coincides with the IR source IRAS 16487- 
4423, which has colours t ypical of ultracompact Hit regions 
(IWood & ChurchwelllfT^ . 

A composite image of S 24 and its environs at 3.6, 4.5, and 
8.0 ;um is displayed in Fig. [T] At 3.6 pm, the emission arises 
in a faint PAH feature at 3.3 pm and in dispersed stellar light, 
while at 8.0 pm most of the emission originates in strong fea¬ 
tures from PAHs. The 4.5 pm band shows emission from Bra 
and Pf/3, and vibrational H 2 lin es and ro-vibrationa l CO lines, 
typical from shocked gas (see I Watson et al.l l2008h . At 8 pm. 
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Fig. 1. Composite image of the IR dust bubble S 24 and its environs showing IRAC images: 3.6 jim is in blue, 4.5 jum is in green, 
and 8.0 jum is in red. The different features in the environment of S 24 are indicated. 


the bubble presents a full filamentary shell of about 28" in ra¬ 
dius, brighter towards the north and east and displays a sharp 
eastern border. Two small and bright regions of extended emis¬ 
sion can be identified at 8 /tm in the image, named G341.220- 
0.213 at RA,Dec.(J2000) = (16'’52'’T8.25X -44°26'54"), of 
about 10" in radius, and G341.217-0.237 at RA,Dec.(J2000) = 
(16''52"*23.74X -44°27'53"), of about 5" in radius, placed to 
the north and east of S 24, respectively. At this wavelength, the 
shape of G341.220-0.213 resembles a horseshoe opened toward 
the eas t. The presence of an in frared dark cloud (IRDC) identi¬ 
fied bv IPeretto & Fulleil (l2009h as SDC341.194-0.221, is evident 
to the west of the bubble where the emission at 8 jim is lacking. 
It is centered at RA,Dec.(J2000) = (16''52'"15.5X -44°28'19.4") 
and extends almost 2' in dec l ination . Two IRDCs were also 
identified by iPeretto & Fulled (l2009t) : one at RA,Dec.(J2000) 
= (16''52“15.4X ^4°30'6.9") (SDC 341.171-0.241), to the 
south of the bubble, and the other one at RA,Dec.(J2000) 
= (16'’52'"35.7X ^4°28'21.8") (IRDC SDC 341.232-0.268). 
Finally, fainter emission can be detected in the whole region, 
encircling S24, G341.22-0.213, and G341.217-0.237. An al¬ 
most circular structure of about 16" in radius, named from 
hereon G341.210-0.252, is also detected in the 8 jum image at 
RA,Dec.(J2000) = (16''52"*26.5X -44°28'50"), to the southeast 
of S 24. The above mentioned features are indicated in Fig. [1] 
The figure also shows that G341.220-0.213, G341.217-0.237, 
and part of the S 24 bubble coincide with sources at 4.5 jum. 

iHenning et al.l (l200lh observed the region at 1.2 mm using 
the SEST telescope with an angular resolution of 23". They de¬ 
tected cold dust emission from G341.220-0.213 and G341.217- 
0.237, and an extended halo linked to S 24, and estimated that 
the total mass in the region amounts to almost 5000 Mq (adopt¬ 
ing a gas-to-dust ratio of 100 and a distance of 3.7 kpc). 


Previous studies toward the region of S 24 have shown 
that star formation is active. Two spots of maser emis¬ 
sion were detect ed. The first spot is a methanol and wa¬ 
ter maser emitter dBavandina et al.ll2012l : lUrauhart et alJ 120131 
IWalsh et al.l 120141) with peak velocities betwe en ~ -50 km 


an d -38 km s * ( Pestalozzi_e^alJ_ 2005Caswell et al.l 


1201 Ol : IVoronkov et al.l 12014 : Walsh et alJ 12014 ). coincident 
with G341.220- 0.213. The second s pot is also a rn e thano l 
maser emitter (IVassilev et al.l l2008t iBavandina et alJ 12012h 
with peak velocities between ~ -70 km s~* and -16 km 
s ' (iLarionov & Val’ttsl l2008t IVoronkov et al.l l2014l) located 
at RA,Dec.(J2000) (16''52'"15X9, -44°28'36'.'7), coincident 

with the IRDC SDC341.194-0.221. Although a number of 
EGO candidates have be en identified in the S21-S24 complex 
(ICvganowski et al.]|2008h . none of them appear projected in the 
close environs of the IR bubble. 


S 24 is detected as a radio continuum source in the Molonglo 
Galactic Plane Sur vey at 843 MHz (MGPSJ165220-442802, 
iMurphv et al.ll200'^ (synthesized beam = 65"x45") with an in¬ 
tegrated flux density S 0.843 = 0.92+0.04 Jy. 

iBronfman et al.l d 19961) performed CS(2-1) line observa¬ 
tions with an angular resolution of 50", using the SEST tele¬ 
scope. They observed the molecular line in one position to¬ 
wards IRAS 16487-4423 (RA,Dec.(J2000) = (16'‘52'"21.24X - 
44°28'275), detecting CS emission with a velocity of —43.4 km 
s ', indicating the presen ce of dense gas towards the IR source. 
iRusseil & CastetsI d2004l) observed the CO(l-O) and CO(2-l) 
lines with an angular resolution of 42" and 21", respectively, 
towards a position located 43 to the southeast of the S 24 bub¬ 
ble, detecting molecular gas having velocities in the range -40 
km s ' to -37 km s“'. All velocities in this paper are referred to 
the LSR. 
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Circular galactic rotation models predict that gas having ve¬ 
locities of -43.4 km s“' lies at near and far kinematical dis- 
tan ces of 3.6 kpc and 12-13 kpc, respectively (see for exam¬ 
ple [Brand|&^i^[l993t). Si milar distances c an be derived from 
maser velocities. Following iHenning et al.l (1200 ll) we adopt a 
distance of 3.7 kpc for the S24 region and its environment. 
Taking into account a velocity dispersion of 6 km s ' for the 
interstellar gas, we adopt a distance uncertainty of 0.5 kpc. 

Here, we present the very hrst molecular line and dust con¬ 
tinuum study towards S 24 and its environs aimed at determining 
the distribution of the molecular gas and the dust linked to the 
bubble and their correlation between them, analyzing their phys¬ 
ical conditions, masses and ambient densities, and investigating 
the status of star formation in the region. The presence of dense 
gas linked to this region makes it particularly interesting to in¬ 
vestigate the physics and kinematics of Hii regions and bubbles 
in dense media, allowing a better understanding of the evolution 
and star formation proce sses around t hese o bjects. Preliminary 
results were published in lCapna et al.l (l2013l) . 


2. Data sets 

2.1. Molecular line observations 

The molecular gas characteristics were investigated by perform¬ 
ing *^CO(3-2) line observations (at 330.58797 GHz) of a region 
of 3:2 X 3:2 during August 2009, and '2C0(2-1), '3CO(2-l), 
and C**0(2-l) line observations of a region of 5' x 5' obtained 
in October 2010, with the APEX telescope using the APEX-1 
and APEX-2 receivers, whose system temperatures are 150 K 
and 300 K, respectively. 

The half-power beam-widths of the telescope are 30" (for 
the (2-1) lines) and 21" (for the (3-2) lines). The data were ac¬ 
quired with a EET spectrometer, consisting of 4096 channels, 
with a total bandwidth of 1000 km s ' and a velocity resolution 
of 0.33 km s“'. The region was observed in the position switch¬ 
ing mode using the OTE technique with a space between dumps 
in the scanning direction of 9" for the (2-1) transition, and in the 
position switching mode in the (3-2) transition. The off-source 
position free of CO emission was located at RA,Dec.(J2000) = 
(13*33'”10.3X-62°2'41"). 

Calibration was performed using Mars and X-TrA sources. 
Pointing was done twice during observations using X-TrA, o- 
Ceti and VY-CMa. The intensity calibration has an uncertainty 
of 10%. 

The spectra were reduced using the Continuum and Line 
Analysis Single-dish Software (CLASS) of the Grenoble Image 
and Line Data Analysis Software (GILDAS) working groupQ. A 
linear baseline fitting was applied to the data. The rms noise of 
the profiles after baseline subtraction and calibration is 0.35 K 
for the (2-1) and (3-2) transitions. The observed line intensities 
are expressed as main-beam brightness temperatures T^b, by di¬ 
viding the antenna temperature Ta by the main-beam efficiency 
7]mb, equal to 0.72 for APEX-1 and 0.82 for APEX-2 (Vassilev 
et al. 2008). The Astronomical Image Processing System (AIPS) 
package and CLASS software were used to perform the analysis. 


' http://www.iram.fr/IRAMFR/PDB/class/class.html 


2.2. Dust continuum data 

To trace cold dust emission we used far-infrared (EIR) images 
from the Herschel Space Observator]^ and from the APEX 
Telescope Large A rea Survey of the Galaxy (ATLASGAL, 
ISchulleretal]l2009tl 


2.2.1. Herschel data 


The archival data correspond to the Hi-GAL ke y program (Hi- 
GAL .Herschel Infrared GALactic plane survey, iMolinari et al.l 
I 2 OIOI OBSIDs; 1342204094 and 1342204095). 


The data hav e been taken in para llel mode with the in¬ 
struments PACS (Po glitsch et alJl^lOl) at 70 and 160jum, and 
SPIRE (lGriffinet al.ir2010h at 250. 350, and 500 jum. The cho¬ 
sen tile covers a held of 2.2 square degrees and is approximately 
centered at (l,b) - {340°,0°). The angular resolutions for the 5 
photometric bands spans from 8" to 3572 for 10 pm to 500 pm. 
Hersc hel Interactive Processing Environment (HIPE vlfl lOtll 
l2010h was used to reduce the data, with reduction scripts from 
standard processing. To combine the two obsids (scan and cross¬ 
scan) of PACS archive data, we used the HIPE implementation 
script of the map merger called Scanmorphos. This script starts 
from Level 1 frames and creates a map from the data. The ob¬ 
tained pixel sizes are 2"/px to blue and 3"/px to red. 


The zero-level of the PACS and SPIRE archive data is un¬ 
known, which may cause problems, for example to derive dust 
temperatures. The zero-l evel offsets of the SP IRE maps were 
corrected as described in [Bernard et al.l (l2010h . A proper zero- 
level was set bearing in mind that Planck includes the same frac¬ 
tion of the sky as Hi-GAL data and taking into account that the 
two common photometric channels of Planck and SPIRE (857 
GHz ~ 350 pm a nd 545 GHz ~ 5 50 pro) are calibrated with re¬ 
spect to LIRAS (iPiat et al.ll2002h . We used the HIPE implemen¬ 
tation script to run destriper gain corrections and to re-calculate 
the absolute offset via cross-calibration with Planck-HEI data 
(HIPE vl2, updated version by L. Conversi 2013). 


To estimate fluxes, a Gaussian smoothing was applied to con¬ 
volve all images at low resolution of 500 pm (a EWHM of 3572). 


Assuming that an extended source has a spectral index = 4 
(2 for the blackbody emission plus 2 for the opacity of the dust), 
to convert from monochromatic intensities of point sources to 
monochromatic extended sources the obtained fluxes of each 
source were multiplied by 0.98755, 0.98741, and 0.96787 for 
250, 350, and 500 pm, respectively. 


2.2.2. ATLASGAL (jata 

This survey was performed at 870 pm (345 GHz) with an angular 
resolution of 19'/2 ( HPBW), using the L arge Apex Bolometer 
Camera (LABOCA) (ISiringo et alJl200^ . The camera operates 
with 295 pixels at the APEX 12-m sub-millimeter telescope. 


^ Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im¬ 
portant participation from NASA 

^ HIPE is a joint development by the Herschel Science Ground 
Segment Consortium, consisting of ESA, the NASA Herschel Science 
Center, and the HIFI, PACS and SPIRE consortia members, see 
http://herschel.esac.esa.int/HerschelPeople.shtml 
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Fig. 2. Overlay of the IRAC emissions at 8 jum (in red) and 
3.6 fim (in green) and the LABOCA emission at 870 jum (in 
blue). Positions of identified YSOs are indicated with squares 
in different colours: Class I (red) and Class II (blue) sources 
from the WISE database. Class I (yellow). Class II (white), and 
Class I/II (green) sources from the Spitzer database, 2MASS 
sources (cyan), and MSX sources (magenta). Dust clumps from 
Contreras et al. (2013) are indicated. 


2.3. Complementary data 

The millimeter and sub-millimeter data were complemented 
with Spitzer images at 3.6, 4.5, and 8.0 ;um from the Galactic 
Legacy Infrared Mid-P lane Survey Extraordinaire (GLIMPSE) 
(iBeniamin et al.l l2003h . and at 24 pm from the MIPS Inner 
Galactic Plane Survey (MIPSGAL) (ICarev et al.ll2005ll . Radio 
continuum images at 843 MHz from the Molonglo Galactic 
Plane Survey (MGPS2, iMumhv et alJl200'^ with a synthesized 
beam of 65"x 45", and at 1.4 GHz fr om the Southern Galactic 
Plane Survey (iHaverkorn et al.lr2006h (beam size = V.l) were 
used. 

To investigate the existence of young stellar objects (YSO) 
candidates projected onto the region we analyzed infrared 
point sources fro m the Midcourse Space Experiment (MSX, 
Price et al. l2001h . the Two Micron All Sky S urvey (2MASS, 
Cutri et al. l2003li . the Spitzer (lEazio et al.ll2004l) , and the Wide- 
held Infrared Survey Explorer (WISE. IWright et al.ll2010l) . 

3. YSOs projected onto the region 

As mentioned in Sect.l, maser emission linked to star-forming 
regions has been detected in the region. Bearing this in mind, we 
decided to search for YSO candidates using the MSX, 2MASS, 
WISE, and Spitzer point source catalogues. To accomplish this 
we analyzed the IR sources in a region of 3' in radius centered 
at RA,Dec.(J2000) = (16''52'"20X -44°28'00"). The results are 
included in Table 1 and their positions are shown in Eig.|2]with 
different symbols. 

To identify YS O candidates in the M SX catalog we fol¬ 
lowed the criteria of iLumsden et alJ d 2002 h . which takes into ac¬ 


count their loci in the (E 2 i/Eg, E 14 /E 12 ) diagram. Two candidates 
were found, classihed as massive young stellar objects (MYSO), 
one of them linked to G341.220-0.213 and the other one to the 
brightest section of the S 24 bubble. 

To search for YSO candidates in the 2MASS catalog, we 
have constructed an [H-K^, J-H] diagram using only sources 
with photometric quality A or B. Six candidates were identified. 
The two upper panels of Eig. [3] show the 2MASS color-color 
(CC) and color-magnitude (CM) diagrams with the YSO candi¬ 
dates represented by red asterisks. 

Regarding the Spitzer database, w e have performed a search 
for candidates adopting the criteria bv lAllen et al.l (l2004t) . which 
allows to discriminate IR sources according to a Class scheme; 
Class I are protostars surrounded by dusty infalling envelopes 
and Class II are objects whose emission is dominated by accre¬ 
tion disks. Considering only sources detected in the four Spitzer- 
IRAC bands, we have investigated their evolutionary stage by 
analyzing their location in a color-color diagram (middle panel 
of Eig. [3ll following t he color criteria determined empirically 
by lAllen et ^ (l2004l) . Class II candidates lie in the region 
0 < [3.6]-[4.5] <0.8 and 0.4 <[5.8]-[8.0] <1.1, and their IR ex¬ 
cess might be produced by accretion disks around the stellar ob¬ 
jects. Sources outside the domain of Class II objects are iden¬ 
tified as Class I candidate objects, and their IR excess might 
be due to circumstellar envelopes around the star. In our study 
we have identified 14 Class I candidates, 5 Class I/II candidates, 
and one Class II. The middle panel of Eig. [3] shows the position 
of the 20 Spitzer sources in the color-color [3.6]-[4.5] vs. [5.8]- 
[8.0] diagram. In this panel, the area delineated by thin lines 
corresponds to the locus of Class I objects, while the area within 
thick dashed lines corresponds to class II objects. 

Einally, with the aim of identifying candidate YSOs using 
photometric data from the Wide-field Infrared Survey Explorer 
survey, which maps the whole sky in four bands ce ntered at 3.4, 
4.6, 12 , and 22 pm, we have used the criteria of iKoenig et al.l 
(l2012h . using [3.4]-[4.6] vs [4.6]-[12] diagrams, as follows: 
after removing contamination arising from background objects 
like galaxies (very red in [4.6]-[12]), broad-line active galac¬ 
tic nuclei (of similar colours as YSOs, but distinctly fainter) 
and resolved PAH emission regions (redder than the majority 
of YSOs), we identified infrared excess sources demanding that 

[3.4] - [4.6] - 0-1 > 0.25 

[4.6] -[12.0]- 0-2 > 1.0 

where [3.4], [4.6], and [12.0] are the WISE bands 1, 2, and 3, re¬ 
spectively, and (Ti = cr([3.4] - [4.6]) and 1 x 2 = cr([4.6] - [12.0]) 
indicate the combined errors of [3.4]-[4.6] and [4.6]-[12.0] col¬ 
ors, added in quadrature. Class I sources are a sub-sample of this 
defined by 

[3.4] -[4.6] >1.0 

[4.6] -[12.0] >2.0 

The rest of the sources are then Class II objects. Identified Class 
I and Class II candidates are shown in the bottom panel of Eig. [3] 
In this panel, the separation between the location of ClassI and 
Class II sources is indicated by the dashed lines. 

As can be seen from Eig. |2l four candidates (one MYSO, 
one 2MASS source, one SPITZER Class I, and one WISE Class 
I) appear projected onto G341.220-0.213. This fact along with 
the presence of maser emission towards the region confirms that 
this is an active star-forming region. Three candidates coincide 
in position with the S 24 bubble, while only a WISE Class II 
candidate is projected onto G341.217-0.237. Most of the candi¬ 
dates appear to the east of the IRDC SDC341.193-0.221 and are 
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Table 1. YSO candidates projected onto S 24 and its environs. 


MSX candidates 

# [hms\ 

[0'"] 

MSX 

name 

8 pm 

Fluxes[Jy] 
12 pm 

14 pm 

21 pm 




Classification 


1 

1652 21.7 

-A4 27 59.8 

G341.2105-00.2325 

6.2092 

14.360 

16.477 

39.985 




MYSO 


2 

16 52 18.5 

-^4 26 55.0 

G341.2182-00.2136 

2.5740 

3.0962 

4.1452 

14.837 




MYSO 


2 MASS candidates: 












# 

a[h m 3'] 

SV ' "] 

2MASS 

Spitz.er 


FluxesTmagl 



Fluxes [mag] 






name 

name 

J 

H 

Ks 

[3.6] 

[4.5] 

[5.8] 

[8.0] 


3 

16 52 19.32 

-44 26 54.2 

16521932-4426541 


14.586 


WSSIM 






4 

16 52 23.17 

-44 30 44.1 

16522316-4430441 

G341.1778-00.2650 

14.653 

HIIBH 


12.658 

12.258 




5 

16 52 23.50 

-44 30 41.3 

16522350-4430412 

G341.1791-00.2651 

14.177 

HBBEH 

13.266 

12.579 

12.373 

12.227 



6 

16 52 26.83 

-44 30 05.6 

16522683-4430055 

G341.1929-00.2666 

15.733 

15.114 

14.377 

13.308 





7 

16 52 30.43 

-44 28 28.0 

16523043-4428279 

G341.2207-00.2576 

14.625 

14.166 

13.781 

13.862 

13.075 




8 

16 52 33.01 

-44 27 57.1 

16523300-4427570 


15.994 

15.233 

14.591 

13.867 

13.231 




Spitzer candidates: 












# 

a[h m 3] 

S[° ' "] 

Spitzer 

2MASS 


Fluxes [magi 



Fluxes [mag] 


Class 




name 

name 

J 

H 

Ks 

[3.6] 

[4.5] 

[5.8] 

[8.0] 


9 

16 52 07.42 

-44 27 05.8 

G341.1950-00.1902 

16520740-4427055 

13.721 

12.690 

12.025 

10.569 

9.877 

9.300 

8.599 

i/ii 

10 

16 52 13.18 

-44 29 49.6 

G341.1708-00.2324 





13.868 

13.610 

10.924 

9.793 

11 

11 

16 52 14.96 

-44 29 22.0 

G341.1800-00.2316 

16521493-4429218 

15.478 

14.088 

<11.675 

11.695 

11.297 

10.701 

10.078 

VII 

12 

16 52 15.10 

-44 28 08.9 

G341.1960-00.2190 





12.845 

10.867 

9.955 

9.422 

I 

13 

16 52 15.44 

-44 27 17.4 

G341.2077-00.2107 





12.808 

11.475 

10.675 

10.211 

I 

14 

16 52 15.98 

-44 26 37.9 

G341.2172-00.2050 





12.690 

11.839 

11.375 

10.424 

VII 

15 

16 52 17.93 

-44 26 53.0 

G341.2176-00.2121 

16521792-4426524 

<17.104 

<13.974 

<11.036 

7.986 

6.215 

5.140 

4.649 

I 

16 

16 52 19.55 

-44 26 40.7 

G341.2233-00.2137 





12.436 

12.049 

9.074 

7.325 

I 

17 

16 52 22.53 

-44 27 17.7 

G341.2210-00.2270 

16522251-4427175 



14.169 

10.627 

9.801 

9.094 

8.498 

I 

18 

16 52 22.98 

-44 25 36.1 

G341.2437-00.2101 





12.595 

12.343 

11.746 

9.833 

I 

19 

16 52 25.16 

-44 27 52.1 

G341.2186-00.2391 





11.203 

9.876 

8.916 

8.330 

1 

20 

16 52 26.03 

-44 28 34.6 

G341.2111-00.2486 





11.410 

10.145 

9.104 

8.676 

I 

21 

16 52 26.71 

-44 28 01.5 

G341.2195-00.2443 





13.473 

11.655 

10.182 

9.131 

I 

22 

16 52 27.59 

-44 30 34.1 

G341.1884-00.2733 





12.480 

11.695 

11.129 

10.675 

lAi 

23 

16 52 27.78 

-44 29 28.5 

G341.2028-00.2622 

16522776-4429282 

<16.378 

14.500 

12.353 

10.098 

9.342 

8.706 

8.303 

VII 

24 

16 52 29.66 

-44 28 46.3 

G341.2154-00.2590 





9.366 

8.211 

7.294 

6.387 

I 

25 

16 52 29.66 

-44 29 27.4 

G341.2066-00.2663 





12.577 

11.575 

10.803 

10.144 

I 

26 

16 52 33.55 

-44 28 27.0 

G341.2269-00.2646 





13.792 

11.341 

9.639 

8.569 

I 

27 

16 52 34.16 

-44 28 36.2 

G341.2261-00.2676 

16523407-4428356 

12.694 

12.351 

12.242 

11.707 

10.458 

9.395 

8.610 

1 

28 

16 52 34.45 

-44 28 35.3 

G341.2268-00.2681 





12.884 

11.260 

10.448 

10.137 

I 

WISE candidates: 












# 

a[h m 3] 

S[° ' "] 

WISE 

2MASS 


Fluxes [magi 



Fluxes [mag] 


Class 




name 

name 

J 

H 

Ks 

[3.4] 

[4.6] 

[12.0] 

[22.0] 


29 

16 52 12.325 

-44 29 50.65 

J165212.32-442950.6 

16521231-4429505 

<17.338 

13.789 

11.393 

10.100 

9.697 

7.192 

4.469 

11 

30 

1652 13.168 

-44 29 37.35 

J165213.16-442937.3 

16521312-4429365 

<15.988 

<14.795 

13.835 

11.981 

11.412 

6.822 

3.941 

II 

31 

16 52 13.541 

-44 29 26.67 

J165213.54-442926.6 

16521357-4429260 

<16.221 

<12.706 

<10.475 

9.009 

8.647 

7.232 

4.362 

II 

32 

16 52 17.106 

-44 27 46.29 

J165217.10-442746.2 





11.963 

10.817 

6.243 

3.134 

I 

33 

1652 17.911 

-44 28 12.87 

J165217.91-442812.8 





10.671 

9.222 

4.283 

0.900 

I 

34 

16 52 17.961 

-44 26 52.95 

J165217.96-442652.9 

16521792-4426524 

<17.104 

<13.974 

<11.036 

8.084 

4.888 

2.607 

-1.190 

I 

35 

1652 18.812 

-44 27 58.92 

J165218.81-442758.9 





7.601 

7.286 

2.830 

-0.197 

II 

36 

16 52 19.544 

-44 28 23.91 

J165219.54-442823.9 

16522015-4428284 

13.624 

11.609 

10.718 

8.584 

8.013 

3.537 

0.262 

II 

37 

16 52 19.890 

-44 28 48.64 

J165219.89-442848.6 

16521984-4428486 

14.585 

<11.779 

<10.221 

8.933 

8.581 

5.674 

1.653 

II 

38 

16 52 21.228 

-44 28 52.70 

J165221.22-442852.6 

16522121-4428522 

<16.713 

<15.318 

13.912 

11.494 

10.541 

5.786 

2.482 

II 

39 

16 52 21.603 

-44 29 04.46 

J165221.60-442904.4 

16522157-4429046 

<15.992 

<14.358 

13.513 

11.138 

10.539 

5.721 

2.730 

II 

40 

16 52 21.886 

-44 26 02.99 

J165221.88-442602.9 

16522182-4426053 

<17.576 

14.702 

12.012 

10.136 

9.590 

8.005 

4.540 

II 

41 

16 52 22.602 

-44 28 47.58 

J165222.60-442847.5 

16522260-4428476 

<17.062 

14.859 

11.133 

8.627 

7.384 

5.155 

1.507 

I 

42 

16 52 22.696 

-44 28 30.59 

J 165222.69-442830.5 





9.250 

8.748 

4.270 

0.807 

II 

43 

16 52 23.588 

-44 27 54.44 

J165223.58-442754.4 

16522366-4427538 

<14.662 

<13.571 

11.702 

7.575 

4.815 

2.094 

-0.978 

I 

44 

16 52 24.638 

-44 29 03.36 

J165224.63-442903.3 

16522467-4429039 

14.852 

12.931 

12.028 

10.809 

10.439 

5.880 

2.396 

II 

45 

16 52 24.659 

-44 29 23.87 

J165224.65-442923.8 

16522482-4429238 

<15.985 

<15.266 

<14.017 

11.731 

11.302 

7.970 

5.320 

II 

46 

16 52 25.457 

-44 28 45.46 

J165225.45-442845.4 





10.425 

9.442 

4.848 

1.575 

II 

47 

16 52 25.605 

-44 27 14.89 

J165225.60-442714.8 





11.632 

10.663 

7.596 

3.503 

II 

48 

16 52 25.668 

-44 28 52.61 

J165225.66-442852.6 

16512576-4428529 

16.343 

<15.243 

<13.988 

10.862 

9.794 

5.340 

2.220 

I 

49 

16 52 25.880 

-44 28 37.18 

J165225.87-442837.1 

16522577-4428393 

<15.257 

<14.320 

14.104 

11.114 

9.575 

5.175 

2.238 

I 

50 

16 52 26.071 

-44 27 26.80 

J 165226.07-442726.8 

16522616-4427259 

<15.232 

<14.120 

13.194 

11.546 

10.794 

6.597 

3.390 

II 

51 

16 52 28.408 

-44 27 23.43 

J165228.40-442723.4 

16522829-4427235 

<17.943 

14.305 

11.736 

10.336 

9.777 

8.372 

5.300 

II 

52 

16 52 29.731 

-44 28 45.25 

J165229.73-442845.2 

16522964-4428459 

<14.407 

<13.557 

<12.079 

9.495 

8.082 

4.686 

1.089 

I 


projected onto emission at 8 /rm, Finally, two Spitzer candidates 
coincide with the IRDC SDC341.194-0.221 and other three can¬ 
didates appear projected onto the borders of G341.210-0.252. 

It is clear from Table 1 that there is little correspondence 
among 2MASS, Spitzer, and WISE candidate YSOs (coinci¬ 
dences: #24 with #52; and #2 with #15 and #34). Regarding 
Spitzer and WISE catalogues, in spite of the similar wave¬ 
lengths of two of the bands and in the different angular res¬ 
olution of the catalogues, an analysis of the coincidence of 
Spitzer and WISE YSOs identified in Table I shows that 
33% of the WISE and Spitzer sources do not have a coun¬ 
terpart in the Spitzer and WISE, respectively, catalogues. 


Additionally, 35% of the WISE and Spitzer sources coincide 
with Spitzer and WISE sources, respectively, which are un¬ 
detected in at least one band and can not be properly classi¬ 
fied. For the rest 17%, there is coincidennce between Spitzer 
and WISE sources, but their counterparts are not classified 
as Class I or II. 


A more detailed analysis will be performed in connection to 
the molecular dust and cold dust in the next sections. 
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Fig. 3. Upper panels: [H-K^, J-H] and [H-K^, K^] diagrams of 
the 2MASS YSO candidates indicated by red crosses. Green as¬ 
terisks correspond to main sequence stars projected onto the S 24 
bubble (see Sect. 7.2). The dashed lines show the reddening vec¬ 
tors corresponding to MO III and B2 V stars in the CC-diagram, 
and to B2 V stars in the CM-diagram. The squares are placed at 
intervals of ten magnitudes of visual extinction. Middle panel: 
Location of YSO candidates from the Spitzer database. Crosses 
represent Class I and II sources. Bottom panel: YSO candidates 
from the WISE database. Class I and II sources are indicated in 
red and green, respectively. 
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Fig. 4. Composite image of the IRAC emission at 4.5 (in green) 
and 8.0 pm (in red), and the Herschel emission at 250 pm (in 
blue). Contours correspond to the emission at 870 pm and goes 
from 0.5 MJy sr"'. 

4. The distribution of the intersteiiar dust 

The distribution of the interstellar dust was investigated using 
data at different wavelengths in the IR. 

Figure |2] shows a superposition of the IRAC 8 pm (in 
red) and 3.6 pm (in green) emissions, and the ATLASGAL 
870 pm emission (in blue). Emission at 870 pm reveals the 
presence of cold dust. This emission approximately encir¬ 
cles the S24 bubble. The dust c lumps identified at this wave¬ 
length by IContreras et al.l (l2013h are indicated. AGAL341.216- 
00.236 and AGAL341.217-00.212 are bright, coincide with 
G341.217-0.237 and G341.220-0.213, respectiv ely, and corre¬ 
late w ith the emission at 1.2 mm detected by iHenning et al.l 
(l200lh . Weaker emission at 870 pm (AGAL SDC341.196- 
00.221) coincides in position with the IRDC SDC341.194- 
0.221 placed to the west of the bubble following its sharp 
eastern border, and is evident due to the absence of 8 pm 
emission in Fig. |2] Additionally, two cold dust clumps are 
present to the south and southeast of S 24, at RA,Dec.(J2000) 
= (16'‘52'"19.6', -44°30'3875) (AGAL SDC34L 172-00.239), 
and RA,Dec.(J2000) = (16*52'”3L7L, -44°29'2671) (AGAL 
SDC34L211-00.271). AGAL SDC34L219-00.259, located at 
RA,Dec.(J2000) = (16''52'"32X -44°28'30"), coincides with a 
section of the IRDC SDC34L232-0.268 and is being analyzed 
by Vasquez et al.(2015, to be submitted). G34L210-0.252 does 
not coincide with emission at 870 pm. 

Figure |4]displays an overlay of the emission at 870 pm (in 
contours) and the emissions at 4.5 (in green) and 8.0 pm (in 
red) from Spitzer, and 250 pm (in blue) from Herschel. Dust 
clumps identified at 870 pm have their counterparts at 250 pm, 
although the emission at 250 pm, which also shows the distri¬ 
bution of cold dust, appears more extended. The brightest dust 
clumps at 250 pm coincide with G34L220-0.213 and G341.217- 
0.237, and are also detected at 350 and 500 pm (the images are 
not shown here). The S 24 bubble appears clearly surrounded 
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to the east and north by SPIRE emission. Herschel emission at 
these wavelengths is also projected onto the IRDC identihed to 
the west of S 24, wh ich can be dehned as a Herschel bright IRDC 
dPeretto et al.ll201(ih . Faint extended emission at 250 jim also co¬ 
incides with the IRDC SDC341.171-0.241, bordering the weak 
halo of 8 jum emission, and near RA,Dec.(J2000) = (16*52'"28 ', 
^4°29'30"), in coincidence with a Spitzer point source identi¬ 
hed at 4.5 /tm. In summary, the correlation between 250 and 870 
/rm is excellent. 

The PACS emission at 70 and 160 jim differs signihcantly 
from the emission at larger wavelengths (Fig. |5]l. The bright¬ 
est source at these wavelengths coincides with G341.220-0.213. 
G341.217-0.237 also emits at 70 and 160 /rm, while a bright 
source is also detected to the east of G341.217-0.237. Faint 
emission correlates with the S 24 bubble, mainly at 70 ^m. This 
emission delineates the sharp western border of the bubble. The 
IRDC SDC 341.194-0.221 can not be identihed either at 70 

or at 160 Emission at 70 jum is also projected onto 
G341.210-0.252, located at 16''52”'26.5X -44°28'50", but it is 
lacking at 160 //m and higher wavelengths. The presence of 
emission at 70 /tm and its absence at 160 jim suggests the ex¬ 
istence of dust at a higher temperature. 

Dust temperatures Tdust for the region using the PACS im¬ 
ages at 70 jim and 160 were estimated by convolving the 
image at 70 jum down to the angular resolution at 160 jum, ap¬ 
plying color corrections, subtracting a background local to the 
S 24 bubble, and assuming that the emission is optically thin. 

The color-temperature map was constructed as the inverse 
function of the ratio map of Herschel 70 pm and 160/tm color- 
an d-background-correc t edmaps, i.e., Td||s^ = f iTY^ (see details 
in iPreibisch et al.ll2012l : [^lendorf et al.ll201'^ . In the optically 
thin thermal dust emission regime f{T) has the following para¬ 
metric form: 

SAv ^ IQ pm) ^ B,{10pm,T) 

’ Sy{v^\60pm) Bv(160jUot,T)\160/ ^ 

where By{v, T) is the Black Body Planck function and j3d cor¬ 
responds to the spectral index of the thermal dust emission. 
The pixel-to-pixel temperature was calculated assuming fid - 
2. This is a typical value adopted for irradiated regions (see 
iPreibisch et alj|20l^ . 

The obtained dust temperature map is shown in Fig. |6l We 
were able to estimate dust temperatures for G341.220-0.213, 
G341.217-0.237, G341.210-0.252, and the S24 bubble. Dust 
temperatures using this method were almost impossible to ob¬ 
tain for the IRDCs due to their low emission at 70 pm. Values for 
G341.220-0.213 are in the range 26-30 K, while for G341.217- 
0.237 temperatures are about 29 K. As regards to the S 24 bub¬ 
ble, the image reveals a clear gradient in dust temperatures (from 
26 to 60 K); the higher the temperature, the closer to the center 
of the bubble, a fact compatible with the presence of ionizing 
sources linked to S 24. Dust temperatures fo r S 24 are compatible 
with values derived bv lWatson et alJ (1201 Oh for the interior of IR 
du st bubbles based on im ages at 24 and 70 pm from MIPS and 
by lAnderson et alJ (1201 2h for Hit regions using Herschel data. 
Dust temperatures for G341.217-0.237 are about 28 K. 

We also estimated Tdust for G341.220-0.213 from the whole 
set of Herschel data and the LABOCA image by fitting a Planck 
function. Fluxes at 70, 160, 250, 350, and 500 pm were calcu¬ 
lated with HIPE using rectangular sky aperture photometry and 
convolving all images down to the angular resolution at 500 pm. 
For each aperture we measured the average surface brightness 
with an uncertainty derived from the standard deviation of the 



35 30 25 16 52 20 15 10 

Right Ascension 


Fig. 5. Uper panel: Overlay of the PACS emission at 70 pm (con¬ 
tours) and the IRAC 8 pm emission (colorscale). Contours are 
0.3 to 0.9 Jy px“\ in steps of 0.2 Jy px“'; 2 to 7 Jy px“', in steps 
of 1 Jy px *; 9, 11,13 to 33 Jy px"', in steps of 4 Jy px”'. Bottom 
panel: Overlay of the PACS emission at 160 pm (contours) and 
the IRAC 8 pm emission (color scale). Contours are 2 to 7 Jy 
px \ in steps of 1 Jy px 10 to 25 Jy px ', in steps of 3 Jy 
px *; and 30 to 65 Jy px ', in steps of 5 Jy px“'. 


surface brightness within each aperture and the flux calibration 
uncertainty. Derived values are 776.0, 1047.2,499.5, 213.3, and 
69.3 Jy at 70, 160, 250, 350, and 500 pm, respectively. Fig. [T] 
shows the spectral energy distribution (SED), which also in¬ 
cludes the flux density at 870 pm obtained by IContreras et al.l 
(1201 3h (14.63+2.47 Jy), and the best fitting (obtained using a 
greybody), which corresponds to Tdust = 28+1 K. Dust tem¬ 
perature for G341.217-0.237 from a spectral energy distribution 
(SED) is difficult to determine. Due to the poor angular resolu¬ 
tion of the Herschel images at 350 and 500 pm, fluxes estimated 
from these data includes part of S 24, what can fake the dust 
temperature. 

Dust masses can be estimated from the expression 
(lHildebrandll983h 


M^dust ~ 


I 870 ' 


IfgVO B870(7’du.st) 


( 2 ) 
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Fig. 6. Dust temperature map derived from the Herschel emis¬ 
sion at 70 and 160 fim. Grayscale goes from 20 to 50 K. Contour 
levels correspond to 20 to 32 K, in steps of 3 K, 35, 40, 50, and 
60 K. Darkest regions indicate highest dust temperatures. 


Spectrum ot 28.540246 beta= 1,5727397 
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Fig. 7. SED for G341.220-0.213 obtained using Herschel and 
ATLASGAL fluxes. 


where S^^o is the flux density at 870 fim, d - 3.7+0.5 
kpc, K’fnn - 1.0 cm^/gr is the dust opacity per unit mass 
(lOssenkopf & Henningi1l9^ . and BgTotFdust) is the Planck 
function for a temperature Tdust- 

Two processes may contribute to the emission at this wave¬ 
length in addition to the thermal emission from cold dust; molec¬ 
ular emission from the CO(3-2) line and free-free emission from 
ionized gas. The continuum emission contribution at 345 GHz 
due to ionized gas was estimated from the radio continuum im¬ 
age at 1.4 GHz considering that the radio emission is thermal 
(see Sect. 5) using the expression ^345 = ( 345 / 1 . 4 )“°' 5 i. 4 . 
Adopting S 1.4 = 1.2 Jy (see Sect. 5), this contribution amounts 
to about 1% of the emission at 870 /rm. To estimate the contri¬ 
bution from the CO(3-2), we took into account the intensity in 
the '^CO(3-2) line. The contribution of this mechanism is less 
than 1%. Thus, the contribution by both mechanisms is within 
calibration uncertainties. 

Adopting Tdust = 28+1 K for G341.220-0.213, and a flux 
^870 = 14.63+2.47 Jy derived by Contreras et al. (2013) for 
AGAL341.217-00.212, a dust mass Mdustx - 12.9+6.3 M© 
can be estimated for G341.220-0.213. Assuming the same 
dust temperature for G341.217-0.237, and a flux 5 870 = 
25.55+4.11 Jy (for AGAL341.216-00.236), we estimate Mdust 2 
= 22.4+11.0 Mq. Finally, for the IRDC to the west of the 
bubble, we adopt temperatures in the range 10-20 K (see 
iPeretto et alJl2010l) . which, along with a flux S 370 = 12.43+2.14 
Jy for AGAL341.196-00.221, results in Mdust-irdc = 17-57 M©. 
Uncertainties in these values are about 50%. Adopting a gas-to- 
dust ratio of 100, gas masses are 1300 M© for G341.220-0.213, 
2300 M© for G341.217-0.237, and 1700-5700 M© for the IRDC. 

5. The ionized gas 

Figure [8] displays an overlay of the radio continuum image at 
843 MHz (in contours) and the emissions at 8 jim and 24 /rm (in 
colorscale). The image at 843 MHz reveals radio emission coin¬ 
cident with the sources, and in particular with the S 24 bubble. 
The emission at 1.4 GHz (not shown here) correlates with that 
at the lower frequency. We derived a flux density 5 1.4 = 1.2 Jy. 
Taking into account a flux density S 0.843 = 0.92 Jy, a spectral in¬ 
dex a - 0.5 (5v oc v“) can be estimated from the emission at both 
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Fig. 8. Overlay of the radio continuum emission at 843 MHz 
(contours), the IRAC emission at 8 fjm (in red), the MIPSGAL 
emission at 24 gm (in blue), and the IRAC emission at 3.6 jum 
(in green). Contours correspond to 10, 30, and 50 mJy beam"', 
and from 100 to 500 mJy beam"' in steps of 100 mJy beam"'. 


frequencies. The cr-value indicates that at least at 843 MHz, the 
radio emission is optically thick, and that the source is thermal 
in nature. 

Figure [8] also shows that S24, as well as G341.220-0.213, 
G341.217-0.237, and G341.210-0.252 emit at 24 gm, indi¬ 
cating the presence of warm dust linked to the sources, and 
consequently, the existence of excitation sources inside them. 
Emission at this wavelength also coincides with the PDR, as re¬ 
vealed by the emission at 8 gm. In fact, the emission at 24 jum 
from Hii regions has two components: emission from very small 
grains out of thermal equilibrium, detected inside the Hii region, 
and ther mal emission from the PDR from grains in thermal equi¬ 
librium d Anderson et al.ll2bl2h . 
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Fig. 10. CO(2-l) (in black), '3cO(2-l) (in red), and C^*^0(2- 
1) (in blue) spectra obtained by averaging the observed spec¬ 
tra in the region of S 24. Intensity is expressed as main-beam 
brightness-temperature. 


Lower limits for the electron density rie and the ionized mass 
Mi can be derived from the image at 843 MHz. For an esti- 
mated flux density of 0.92 Jv , using the classical expressions by 
(iMezger & Hendersonlll96'7l) . a typical electron temperature of 
10^ K, and bearing in mind a He abundance of 10% by particle 
number and that He is single ionized, we derived > 200 cm“^, 
and Mi > 22 Mq. We note that similar results can be obtained 
using the emission at 1.4 GHz. Clearly, the fact that the emission 
at 1.4 GHz is optically thick also can not be ruled out. Bearing 
in mind that the emission at 8 /im due to PAH in the photodis¬ 
sociation region marks the boundary of the ionized region, we 
measured a radius of 0.33 pc for the Hit region. 

Although these results can not be used to determine either the 
ionized mass or the electron density, they reveal the existence of 
an Hii region. 

6. The distribution of the moiecuiar gas 

Figure |9] displays the CO(2-l) and '^CO(2-l) spectra obtained 
for the region under study. The spatial separation between these 
prohles is 30". Relative coordinates are expressed in arcseconds, 
referred to RA,Dec.(J2000) = (16''52'"20X -44°27'30"). The in¬ 
dividual spectra show that the molecular emission appears con¬ 
centrated between —55 km s“' and -20 km s“', where a number 
of velocity components are present, and an isolated component 
near -10 km s ' (which is not going to be taken into account). 
The CO(2-l)/'^CO(2-l) line ratios are close to 1-1.5 in the cen¬ 
tral and northeastern parts of the hgure, suggesting that the gas 
is optically thick at '^€0(2-1) where S 24, G341.220-0.213, and 
G341.217-0.237 are located. 

The CO spectra of Fig. [10] illustrate the general behaviour 
of the molecular gas towards the entire region. There, we show 
the CO(2-l), '^CO(2-l), and C'^0(2-l) spectra averaged within 
the observed region. At least three velocity components can be 
identihed in the line of sight within the velocity interval from 
-55 km s“' to -20 km s“'. 

The bulk of the molecular emission is detected in the three 
molecular lines and peaks at ^-45 km s', while fainter gas com¬ 
ponents are present at ^-38 km s ' and -30 km s '. The faintest 
component is present near-10 km s ' in the '^CO(2-l) line only. 
The velocity components peaking at -38 km s ' and -30 km s ' 


can be barely identihed in the '^CO(2-l) spectrum and are ab¬ 
sent in the C'^0(2-l) prohle. Gas components having the same 
velocity can be also identihed in the averaged J = 3-2 prohle (not 
shown here). _ 

Circular galactic rotation models (e.g. iBrand & BlitzIflOQ^ 
predict near and far kinematical distances in the range 3.2-3.9 
kpc and 12-13 kpc, respectively, for the gas component at -45 
km s ', while near and far kinematical distances in the range 
2.4-3.5 and 12-14 kpc, respectively, were obtained for gas in the 
components at -38 km s ' and -30 km s '. Particularly for the 
component at ^5 km s ', near kinematical distances are com¬ 
patible with previously adopted values for the region of S 24. 

To illustrate the spatial distribution of the molecular gas we 
show in Fig. [TT] the emission of the '^CO(2-l) in the inter¬ 
val ^9.9 km s“' to -40.5 km s“', in steps of 1.05 km s“' (in 
grayscale) superimpossed to the emission at 8 jum (in contours) 
for comparison. Most of the molecular emission has velocities 
in the interval from ^6.8 km s“' to ^1.5 km s"'. 

As regards the S 24 bubble, CO emission encircles the east¬ 
ern and northern sections at -46.8 km s ' and -45.7 km s '. In 
the images at ^4.6 km s“', -43.6 km s“', and ^2.5 km s“' this 
emission extends to the western section of the bubble, where the 
IRDC SDS 341.194-0.221 is present. S 24 appears almost com¬ 
pletely encircled by molecular gas in the images at ^3.6 km s ' 
and -42.5 km s“', while at ^1.5 km s“' and ^0.5 km s“' the 
CO emission is present to the west of the bubble, coincident with 
the IRDC. 

A close inspection of the emission in the interval from ^8.9 
km s“' to ^2.5 km s“' reveals bright molecular emission linked 
to G341.217-0.237 and the eastern border of the IR bubble. 

The CO emission associated with the different IR sources 
described above will be analyzed in more detail in the following 
sections. Material linked to the IRDC SDC341.232-0.268 and 
present in the interval from ^9.5 km s“' to ^2.1 km s“' is 
being analyzed in detail by Vasquez et al. (in preparation). 

7. Putting all together 

In this section we investigate in more detail the gas distribution 
in S 24 and the different features in its environment and its rela¬ 
tion to the warm and cold dust and the YSO candidates. 

7.1. G341.220-0.213 

Panel A of Fig. [T2| displays the '^CO(3-2) emission having ve¬ 
locities in the range ^3.8 km s“' and -42.0 km s“' superim¬ 
posed onto the emission at 8 fim for comparison. The emission 
distribution at these velocities reveals a molecular arc-like struc¬ 
ture closely encircling the bright northern and western borders 
of G341.220-0.213. At 8 fim, G341.220-0.213, which resembles 
a horseshoe, is open opposite to the molecular emission. The IR 
source is clearly interacting with the surrounding molecular gas 
and is evolving in a medium with a density gradient. The coin¬ 
cidence of G341.220-0.213 with emission at 24 yum (see Sect. 5) 
suggests the presence of excitation sources inside. Bearing in 
mind the morphology of the source at 8 um an d in molecular 
lines, a champagne flow dGarav & Lizanol 1999h can not be dis¬ 
carded. 

A number of YSO candidates appears inside the horseshoe: 
one MSX source classified as MYSO (# 2 in Table 1); one 
Spitzer and one WISE sources coincident in position, both clas- 
sihed as Class I (# 15 and # 34). We inspected the spectral energy 
distribution (SED) of # 2, # 15, and # 34, assumed to be only 
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Fig.9. Decl. vs. R.A. map showing *^CO(2-l) (in red) and '^CO(2-l) (in black) profiles observed in a 5'x5' region centered in the 
S 24 bubble. Each profile shows Tmb in the interval from -1 K to 20 K vs. LSR velocity in the range from -80 km s ' to 0 km s '. 
The (0,0) position corresponds to RA,Dec.(J2000) = (16''52"‘20^, ^4°27'30"). Note that R.A. inceases towards the right. 



Fig. 15. Overlay of the '^CO(3-2) emission distribution at - 
47.4 km s"' (in contours) and the emission at 8 /tm (in 
grayscale) showing the molecular clump at the interface between 
G341.217-0.237 and the S 24 bubble. Contours are from 3 K to 
6 K (Tnib), in steps of 1 K. The symbols are the same as in Fig. 
12 . 


one source, using the on-lin^l tool developed bv iRobitaille et al.l 
http://caravan.astro.wisc.edu/protostars/ 


(l2007h . which can help to discriminate between evolved stars 
and reliable candidate YSOs. This tool fits radiation transfer 
models to observational data according to a minimization al¬ 
gorithm. Models that accomplished the following condition were 
selected; 

- aL < 3n (3) 

where is the minimum value of the x^ among all models, 
and n is the number of input data fluxes. 

To perform the fitting we used the photometric data listed 
in Table 1, along with fluxes derived from Herschel and 
ATLASGAL images, visual absorption values in the range 30- 
40 mag, and distances in the range from 3.2 to 4.2 kpc. The re¬ 
sult is shown in the left panel of Fig. [T3] Fluxes obtained from 
Herschel and ATLASGAL were considered upper limits. The 
best fitting suggests a massive central source of about 10 Mq, 
a disk mass of 5x10“"' Mq, an envelope mass of 21 Mq, a total 
luminosity of 5300 Lq, and an age of 1x10® yr. With this char- 
ac teristics, the evo l utiona ry stage of the central source according 
to lRobitaille et akl (l2007h would correspond to Stage II (objects 
with optically thick disks). These results points to a massive cen¬ 
tral object. 

We propose that this source has started modifying its envi¬ 
rons by dissociating and then ionizing the nearby material, orig¬ 
inating a PDR detected by its emission at 8 jim at the interface 
between the ionized and molecular material. It may have con¬ 
tributed to the heating of the dust as revealed by the emission at 
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Fig. 11. Channel maps of the '^CO(2-l) line (in grayscale) for the velocity interval from -4-9.9 km s ^ to ^0.5 km s * in steps of 
1.05 km s * overlaid onto the 8 fim emission (in contours). Grayscale goes from 0.5 K to 17.0 K. 


24 jum. Also, the detection of maser emission towards G341.220- 
0.213 (see Sect. 1) indicates that star formation is active. 

A low emission region is present at RA,Dec.(J2000) = 
(16''52"'22.5X -44°26'40"), probably linked to a weak IR arc 
seen slightly to the north of that position. 

Panels B and C of Fig. [12] show the emission distribution in 
the ^^CO(2-l) and C'®0(2-1) lines for the same velocity interval 
as panel A, while in panel D, the emission corresponds to the 
CO(2-l) line in a slightly different velocity interval. The emis¬ 
sion distribution reveals a molecular clump coincident with the 
G341.220-0.213. No signs of the arc-like structure are present 
in the (2-1) transitions, probably due to the lower angular res¬ 
olution of these lines (30") in comparison with the one of the 
13C0(3-2) (20"). 


The molecular clump detected in the interval [^4.0,^2.0] 
km s“' coincides with the emission of cold dust as imaged in the 
far-lR by LABOCA (AGAL 341.217-00.212) and Herschel (see 
Fig ©, indicating that gas and dust are well mixed in the region. 

The analysis of the data cubes allows us to identify molecular 
gas associated with G341.220-0.213 in the velocity interval from 
^6.5 km s“* to ^0.2 km s“'. 


7.2. The S24 Hii region 

The upper left panel of Fig. [T4l shows the '^CO(3-2) emission 
distribution at ^3.2 km s“* in grayscale and contours, while the 
upper right panel displays an overlay of the same CO contours 
and the emission at 8 ;um. The middle and bottom panels show 
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Fig. 12. Panels A-D: Overlay of the CO emission distribution in contours and the emission at 8 /rm in grayscale towards G341.220- 
0.213. Panel A: '^CO(3-2) line emission in the interval from ^3.8 km s“' to ^2.0 km s“'. Contours are from 5 K to 12 K (Tnib) in 
steps of 1 K. the different symbols mark the position of YSO candidates: stars: MSX sources, triangles: 2MASS sources, crosses: 
Spitzer sources, diamonds: WISE sources. Panel B: '^CO(2-l) line emission in the interval from ^4.0 km s ' to ^1.9 km s“*. 
Contours are from 5 K to 11 K (Emb) in steps of 1 K. Panel C: C*^0(2-l) line emission in the interval from ^3.8 km s * to ^1.9 
km s“'. Contours are from 1.5 K to 4 K (Emb) in steps of 0.5 K. Panel D: *^CO(2-l) line emission in the interval from ^5.7 km s"' 
to ^4.3 km s“'. Contours are from 9 K to 14 K (Emb) in steps of 1 K. 



0.1 1 10 100 1000 0.1 1 10 100 1000 
A (/im) A {fim) 


Fig. 13. Left panel: SED of YSOs from Table 1 in G341.220-0.213. Input data are indicated by filled circles. The triangles indicate 
upper limits. The black line corresponds to the best fitting. The fittings which obey Eq. 1 appear indicated by gray lines. The dashed 
line shows the emission of the stellar photosphere including foreground interstellar extinction. Right panel: SED of source #43 from 
Table 1 in G341.217-0.237. 
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Fig. 14. Upper panels: Left: '^CO(3-2) emission distribution at ^3.2 km s“' showing the molecular envelope around the S24 
bubble. Contours are 3.0 K to 6.0 K (Tnib), in steps of 1.0 K; 8.0 K, 10.0 K, and 12.0 K. Right: Overlay of the same contours of 
the left panel and the IRAC emission at 8 pm. The symbols have the same meaning as in Fig. 12. Middle panels: Left: '^CO(2-l) 
emission distribution at ^3.2 km s ^ Contours are from 4.0 K to 12.0 K in steps of 1.0 K. Right: Overlay of the same 

contours of the left panel and the IRAC emission at 8 pm. Bottom panels: Left: C'^0(2-l) emission distribution at ^3.2 km s '. 
Contours are from 1.0 K to 4.0 K in steps of 0.5 K, and 5 K. Right: Overlay of the same contours of the left panel and the 
IRAC emission at 8 pm. 


the emission at the same velocity for the '^CO(2-l) and C'^0(2- 
1 ) lines, respectively. 

The emission at this velocity reveals an U-shaped, al¬ 
most complete molecular shell encircling the S24 bubble. The 
strongest emission regions of this molecular shell encircles the 
brightest sections of S 24 as seen at 8 pm. The shell appears open 
towards RA,Dec.(J2000) = (16''52“19", -44°28'30"). The emis¬ 
sion in the western part of the shell coincides with the IRDC 
SDC341.194-0.221, which borders this section of the bubble. 
The molecular shell is present in the interval from -A6 km s ' 
to ^1.5 km s“'. 


As shown in FigHT] in the velocity interval from ^6.5 km 
s ' to ^3.5 km s ' the molecular gas borders the eastern, north¬ 
ern and western sections of S 24. Molecular gas at these ve¬ 
locities coincides with cold dust detected in the far IR at d > 
160 pm, as shown in Figs. 2 and 4. Molecular shells have been 
found around many IR dust bubbles, in particular those identihed 
as Hn regions (e.g. iPetriella et al.ll2010l : iPomares et al.l[^09bt 
iDeharveng et al Jl201ol) . The presence of CO in the Hii region in¬ 
terior suggests that the shell corresponds to a three-dimensional 
structure. The appearance of the shell indicates that it is unho- 
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mog eneous, with the mol ecular gas distributed unisotropic ally 
te.g. lAnderson et ^l2015h . 

The molecular shell might be expanding at 2-3 km s ' due 
to the difference in pressure between the ionized gas inside the 
bubble and the neutral gas outside, although no signs of expan¬ 
sion are detected in position-velocity diagrams. We note, how¬ 
ever, that low expansion velocities as those expected in this case 
are not easy to detect. In the case of expansion, the emission de¬ 
tected towards the interior of the region would correspond to the 
approaching and receding sections of the shell. 

Three YSO candidates are projected onto the S 24 bubble: 
one MYSO (# 1 in Table 1), and two WISE sources identified as 
Class II (# 35 and # 36). 

Also, a number of YSOs coincide in position with the molec¬ 
ular shell, suggesting that star formation is active in the environs 
of the bubble. 

Two YSOs appear projected onto the IRDC SDC341.194- 
0.221. They do not coincide with the spot of maser emission 
described in Sect. 1 projected onto this IRDC. 

In addition to these YSO candidates, three stars could be 
identified in the 2MASS catalog in the region under study. 
These stars are indicated with green asterisks in the [H- 
Kj, J-H] and [H-Ks, Ks ] diagrams of Fig. [3] In particular, 
2MASS J17502070-4428012 is projected inside the bubble at 
RA,Dec.(J2000) = (16*52'”20.7X -44°28'T/2). It locus in the 
CC- and CM diagrams (H-K,) = 1.19 mag, (J-H) = 2.18 mag, 
and Ks = 10.94 mag) is indicative of a massive star with an opti¬ 
cal absorption Ay - 30 mag. Although more studies are necces- 
sary, this star might be responsible for dissociating and ionizing 
the gas, creating the observed PDR at 8 jum. 

7.3. The interface between S24 and G341.217-0.237 


are heating the dust and probably dissociating and ionizating the 
molecular gas. 


8. The parameters of the molecular gas 

The molecular mass linked to G341.220-0.213, the clump at the 
interface between G341.217-0.237 and S 24, and the S 24 Hii re¬ 
gion, and the molecular mass in the whole region were estimated 
from the emission in the J=2—> 1 line. The derived parameters are 
listed in Table|2] 

Assuming Local Thermodinamic Equilibrium (LTE) condi¬ 
tions and that the emission in the CO(2-l) line is optically thick, 
we derived the excitation temperature T^xc from the emission in 
the CO(2-l) line using 




ni \ 



-1- 

Ep-izco - 7’i2co 


gl-exc - 1 

- 

e^H - 1 



where T*^ = hv^/k, vn the frequency of the '^CO(2-l) hue, 
and Tbg = 2.7 K. To obtain the peak main beam tempera¬ 
ture in the CO(2-l) hne (Tp nco) we took into account the 
averaged spectrum within regions of about 30" in diameter 
for the cases of G341.220-0.213 and the interface between 
G341.217-0.237 and S 24, while for the S 24 region and the 
whole region the spectra were averaged in a region of 2' and 
3' in diameter, respectively. 

The optical depth T 13 was obtained from the *^CO(2-l) line 
by assuming that the excitation temperature is the same for 
CO(2-l) and '^CO(2-l) emission lines using the expression 


T13 


-In 
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(5) 


As described in Sect. 6 , a bright molecular clump appears pro¬ 
jected onto the interface between G341.217-0.237 and S24 at 
RA,Dec.(J2000) = (16'’52'"23X -44°28'). Figure [B] displays an 
enlargement of the region of the clump showing the emission in 
the *^CO(3-2) line at ^7.4 km s“* (in contours) and the 8 ;um 
image in grayscale. The peak CO emission is coincident with the 
interface between the two sources. The anticorrelation between 
the CO and IR emissions strongly suggests that the molecular 
clump has modeled the eastern border of the PDR in the Hii re¬ 
gion. 


These facts suggest that S 24 and G341.217-0.237 are inter¬ 
acting with the molecular clump. The molecular clump coincides 
with emission at 870 fim and in the Herschel bands at 250 and 
350 jum. 

The YSO candidate # 43, classified as Class I, coincides with 
G341.217-0.237. The SED corresponding to this source is shown 
in the right panel of Fig. [T3] To perform the fitting we used the 
photometric data listed in Table 1, and the flux derived from the 
image at 70 fj.m. Herschel fluxes at larger wavelengths include 
a significant part of the S 24 bubble and were not used for the 
SED. Visual absorption is in the range 30-50 mag, and distances 
from 3.2 to 4.2 kpc. The best fitting corresponds to a central 
source of about 15 M©, a disk mass of 8x10“^ M©, an envelope 
mass of 1.1x10“^ M©, a total luminosity of 22000 L©, and an 
age of 1.3x10® yr. According to this fitting, the object is prob¬ 
ably in Stage III (objects with optically thin disks). However , 
as pointed out bviDeharveng et akl (l2012bll . [Offner et^ (l2012h . 
and iRobitaillel (12008 ), these results should be taken with caution. 
The presence of emission at 24 /urn coincident with G341.217- 
0.237 suggests that the ionizing photons from the massive source 


where 7*^ = hvi^lk, V 13 the frequency of the *^CO(2-l) line. 
Assuming LTE, the *^CO column density, NtoiC^CO), can be 
estimated from the *^CO(2-l) line data following 


NtotC^CO) = 3.2x10*'' 




1 — e 



dv (cm- 2 )( 6 ) 


The integral of Eq.j^can be aproximated by 




(7) 


This aproximation helps to eliminate to so me extend optical 
depth effects and is good within 15% for r < 2 (iRohlfs & WilsonI 
I 2 OO 4 I 1 . Bearing in mind the Ti 3 -values listed in Table 2, it is ap¬ 
propriate for our region. Then, the molecular mass was calcu¬ 
lated using 

M(H2) = (mx,„)-^ fi A N(H2) d^ (M©) (8) 

where msun is the solar mass (~2xl0^^ g), ^ is the mean 
molecular weight, which is assumed to be equal to 2.76 af¬ 
ter al lowance of a relative helium abundance of 25% by 
mass (lYamaguchi et aljll999l) . mn is the hydrogen atom mass 
(~ 1.67x10“^'' g), A is the solid angle of the CO emission, d is the 
adopted distance expressed in cm, and A(H 2 ) is the H 2 column 
density, obtaine d using an abundance A(H 2 ) / V(*^CO) = 5x10® 
(lDickmanll978h . Uncertainties in molecular masses and ambient 
densities are about 50% and 70%, respectively, and originates 
mainly in distance uncertainties and optically thick emission in 
the '®CO(2-l)line. 
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Fig. 16. LVG plots for G341.220-0.213, G341.217-0.237, and the IRDC SDC341.194-0.221. Green lines mark the solutions for 
the ratios '^CO(3-2)/'^CO(2-l), while red lines, for the ratios '^CO(2-l)/CO(2-l). Central lines correspond to the observed ratios. 
Outer lines indicate ~10% calibrations uncertainties. 


Table 2. Molecular gas parameters 



I("CO) 

[K km s-*] 

Texc 

ri3co 

[K] 

R 

["= pc] 

A("CO) 

[10‘* cm-2] 

MHz) 
[10** cm-*] 

^mol 

[Mo] 

m 2 

[10* cm-*] 

G341.220-0.213 

52.3 

19.1 

1.95 

27 = 0.44 

14.6 

7.3 

1550 

6.6 

Interface G341.217-0.237/S24 

18.9 

21.0 

1.65 

10 = 0.16 

5.2 

2.6 

470 

41.6 

S 24-shell 

34.4 

16.8 

1.12 

48 = 0.77 

7.3 

3.6 

3100 

2.5 

Whole region 

32.4 

16.1 

1.14 

114= 1.8 

7.0 

3.5 

10300 

0.59 


To explore more robustly the issues of temperature and 
density on G341.220-0.213 and G341.217 -0.237 we performed 
a large velocity gradient (LVG ) analysis IScoville & SolomonI 
Il973t iGoldreich & Kwanl[l974h for radiative transfer of molec¬ 
ular emission lines. We use the Ivg task implemented as part 
of the MIRIAE0 package of SMA. For a given column den¬ 
sity (normalized by the line width), this program estimates the 
line radiation temperature of a molecular transition as a func¬ 
tion of the kinetic temperature (Tk) and the H 2 volume den¬ 
sity {tiHi). We generated 50 x 50 model grids for the '^CO(3- 
2), '^CO(2-l), and CO(2-l) lines with Tk in the range 1-300 
K and nni with the range 1x10^ to 1x10* cm“*. In all cases 
we convolved the '*CO(3-2) data with a Gaussian of FWHM 
30" in order to smooth the cubes down to the angular resolu¬ 
tion the '*CO(2-l) data. We then derived peak main beam tem¬ 
peratures and line widths for the studied components by fitting 


^ http://www.cfa.harvard.edu/sma/miriad/packages/ 


Gaussian functions to the line profiles. Two ratios, R( 3 _ 2 )/( 2 -i) = 
r„b‘^CO(3-2)/CO(2-l) and Ri2/i3=r„b‘^CO(2-l)/CO(2-l), are 
then calculated. The input column densities of CO and '*CO 
were estimated from the dust-derived column density measured 
from the 870 pm emission in each region, and assuming rela¬ 
tive abu ndances N(CO)/N (H 2 ) = IxlO^^ and N('*CO)/N(H 2 ) = 
2xl0~® (lDickmanlll978l) . The same procedure was applied to 
the IRDC SDC341.194-0.221, for which we have not derived 
molecular mass and ambient density since its molecular emis¬ 
sion is included in the S 24 molecular shell. 

FigHS] shows the solutions for G341.220-0.213, G341.217- 
0.237, and the IRDC. They indicate that kinetic temperatures 
in the I RDC are about 2 0 K, not different than derived in other 
IRDCs dEgan et al.ll998h . and not much higher than typical tem¬ 
peratures of molecular gas without an extra heat source other 
than interstellar radiation field. The volume density of the IRDC 
is about 10^ cm-*. For the case of G341.220-0.213, G341.217- 
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0.237, their temperatures are a bit higher (~ 30 - 40 K), not 
surprising given that they are active star-forming regions (see 
Sect. 7). Regarding their volume densities, they are less well 
constrained but values higher than 4x10"^ cm“^ are inferred from 
the plots. 

An inspection of Table 2 shows that molecular gas densi¬ 
ties are in the range 10"^-10^ for particular regions including the 
shell around the S 24 Hit region. Masses derived from molecu¬ 
lar data are compatible, within errors, with those obtained from 
the emission at 870 ^m, assuming a gas-to-dust ratio equal to 
100. Additionally, molecular ambient densities are compatible 
with values derived from LVG analysis. The large difference for 
G341.217-0.237 arises in the fact that the mass estimate from 
molecular emission corresponds only to the emission at the inter¬ 
face between G341.217-0.237 and the S 24 Hit region. The am¬ 
bient density obtained by distributing the total molecular mass 
in the whole observed region is high (5900 cm“^), indicating 
that G341.220-0.213, G341.217-0.237 and S 24 are evolving in 
a high density ambient medium. 

The high opacity estimated for G341.220-0.213 is not sur¬ 
prising bearing in mind the high column density of this source. 


9. Triggered star formation scennario in S24? 

We have found evidences for ongoing star formation in the 
dense molecular envelope encircling the S 24 Hit region. Have 
these YSOs been triggered through the collect and collapse (CC) 
mechanis m? To answer this questi on we can apply the analytical 
model by IWhitworth et al.l (Il994h . For the case of Hit regions, 
the model predicts the age of the Hit region at which the frag¬ 
mentation occurs (the fragmentation time scale), f^ag, the size of 
the Hit region at that moment, Rftag, the mass of the fragments, 
Mfrag, and their separation along the compressed layer, r^ag. The 
parameters required to derive these quantities are the UV pho¬ 
ton flux of the exciting star, ALy, the ambient density of the sur¬ 
rounding medium into which the Hit region is evolving, no, and 
the isothermal sound speed in the shocked gas, a^. 

Since the exciting star is unknown, we took into account a 
large range of spectral types, i.e. from 05V to 0 9V stars, with 
UV fluxes in the range ALy = (18-0.8)xl0"^* s“' (iMartins et alJ 
l2002h . Using the mean H 2 ambient density nni - 5900 cm“^ 
(see Table 2), and a^ = 0.2-0.6 km s“', we obtained ffrag = (6.6- 
8.8)xl0^ yr, Rg-ag = 1.8-2.3 pc, M^ag = 22-29 Mq, and, r^ag = 
0.5-0.4 pc. 

The dyna mical age of the Hn regio n can be estimated using 
the equation (iDvson & William^ll997h 


Uyn 


7cs Us/ 



(9) 


where Rs is the original Stromgren radius, equal to 0.09-0.25 
pc for the adopted spectral types, and Cs is the sound velocity 
in the ionized gas. Derived dynamical ages span the range (9- 
43)xl0^ yr. We hnd that the dynamical age is smaller than the 
fragmentation time scale f^ag for the adopted ambient density. 
The comparison of these two quantities does not support a col¬ 
lect and collapse process, which is a relatively slow process, for 
the triggering of star formation in the envelope. An RDI sce¬ 
nario could be investigated, however evidences of this process 
(such the presence of pillars) appear to be absent. Very probably, 
the compact Hit region is too young for triggering to have begun. 


10. Conclusions 

Based on the molecular emission in the CO(2-l), ^^CO(2-l), 
C**0(2-l), and ^^CO(3-2) lines, obtained with the APEX tele¬ 
scope, and images in the near-, mid-, and far-IR from IRAC- 
Glimpse, Herschel, and ATLASGAL we performed a multi¬ 
wavelength analysis of the infrared dust bubble S 24, and two 
extended IR sources (G341.220-0.213 and G341.217-0.237) lo¬ 
cated in its environs. The region coincides with the IRAS point 
source 16487-4423, classihed as UCHII. We also investigated 
the presence of YSO candidates in the region using the MSX, 
2MASS, Spitzer, and WISE catalogs. 

The molecular emission distribution shows that gas linked 
to the S24 bubble, G341.220-0.213, and G341.217-0.237 has 
velocities between ^8.0 km s“' and ^0.0 km s“\ compatible 
with a kinematical distance of 3.7 kpc previously adopted for the 
region. 

The gas distribution reveals a shell-like molecular structure 
of ~0.8 pc in radius detected in the velocity interval from -Ab 
km s“' to ^1.5 km s“', encircling the S 24 bubble. A cold dust 
counterpart of this shell is detected at wavelengths larger than 
160 pm, i.e. in the LABOCA and SPIRE images. On the con¬ 
trary, PACS emission at 70 pm, as well as MIPS emision at 24 
pm from warm dust appear projected inside the bubble. Dust 
temperatures derived from the emission at 70 and 160 pm are 
in the range 26 to 60 K, whith the higher values closer to the 
center of the bubble. The detection of radio continuum emis¬ 
sion and the presence of warm dust indicates the existence of 
exciting sources and that the bubble is a compact Hit region. We 
estimated a molecular gas and H 2 ambient density in the shell 
of 3100 Mq and 2.5x10"^ cm“^. A search for excitation sources 
allows the identihcation of a massive star candidate projected 
onto the Hit region, which would provide the ultraviolet photons 
to heat the gas and dissociate and ionize the gas. A number of 
YSO candidates appear projected onto the molecular shell. 

Part of the molecular gas bordering the S 24 bubble coincides 
with the extended IRDC SDC341.194-0.221, for which we de¬ 
rived a kinetic temperature of about 20 K and a volume density 
of 10"^ cm“^, based on LVG analysis. 

As regards G341.220-0.213, which resembles a horseshoe 
open towards the east, the '^CO(3-2) line emission allowed to 
detect an arc-like structure with velocities in the range ^3.8 
km s“* to ^2.0 km s“*, encircling the brightest sections of 
the IR source. The morphology indicates that the molecular 
gas is interacting with the IR source. A dust counterpart is de¬ 
tected in the Herschel (SPIRE and PACS) and LABOCA images. 
The spectral energy distribution constructed using Herschel and 
ATLASGAL fluxes allows to derive a dust temperature of 28+1 
K. A molecular mass of 1550 Mq is linked to this source, com¬ 
patible with the mass derived from ATLASGAL. The H 2 am¬ 
bient density amounts to 6.6x10"^ cm“^, in agreement with the 
ambient density value estimated from LVG analysis. A number 
of IR point sources having characteristics of YSO coincide with 
this region, which along with the presence of maser emission 
indicates that star formation is active. The SED for the derived 
fluxes suggests the presence of a massive central source of 10 
Mq, which would be responsible for heating the dust and disso¬ 
ciating the molecular gas. 

G34L217-0.237 also coincides with molecular gas and cold 
dust. The line emission allowed the detection of a molecular 
clump present at the interface between the S 24 Hit region and 
G341.217-0.237, shapping the eastern border of the IR bubble. 
A high ambient density (4.2x10^ cm“^) was estimated for the 
clump. 
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The total mass of molecular gas in the region and the H 2 
ambient density amount to 10300 Mq and 5900 cm“^, indicating 
that G341.220-0.213, G341.217-0.237, and the S24 Hit region 
are evolving in a high density region. 

Finally, we investigated if a collect and collapse scenario can 
explain the presence of YSOs in the dense envelope around the 
S 24 Hit region. We conclude that neither this mechanism nor the 
RDI process seem to occur in this region. Very probably, the Hit 
region is too young for triggering to have begun. 
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